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ABSTRACT 

We study the Sloan Lens ACS survey (SLACS) strong lensing system SDSS 
J1430+4105 at Z\ = 0.285. The lensed source (z s — 0.575) of this system has a complex 
morphology with several subcomponents. Its subcomponents span a radial range from 
4kpc to lOkpc in the plane of the lens. Therefore we can constrain the slope of the 
total projected mass profile around the Einstein radius from lensing alone. We mea- 
sure a density profile that is slightly but not significantly shallower than isothermal 
at the Einstein radius. We decompose the mass of the lensing galaxy into a de Vau- 
couleurs component to trace the stars and an additional dark component. The spread 
of multiple image components over a large radial range also allows us to determine 
the amplitude of the de Vaucouleurs and dark matter components separately. We get 
a mass to light ratio of Md £^ auc w (5.5 ± 1.5) and a dark matter fraction within 

the Einstein radius of w 20% to 40%. Modelling the star formation history assuming 
composite stellar populations at solar metallicity to the galaxy's photometry yields a 
mass to light ratio of « 4.01? j ^ and « 2.3t°j ^ for Salpeter 

and Chabrier IMFs, respectively Hence, the mass to light ratio derived from lensing 
is more Salpeter-like, in agreement with results for massive Coma galaxies and other 
nearby massive early type galaxies. We examine the consequences of the galaxy group 
in which the lensing galaxy is embedded, showing that it has little influence on the 
mass to light ratio obtained for the de Vaucouleurs component of the lensing galaxy. 
Finally, we decompose the projected, azimuthally averaged 2D density distribution of 
the de Vaucouleurs and dark matter component of the lensing signal into spherically 
averaged 3D density profiles. We can show that the 3D dark and luminous matter 
density within the Einstein radius (R,Ein ~ 0.6R o ff ) of this SLACS galaxy is similar to 
the values of Coma galaxies with the same velocity dispersions. 

Key words: gravitational lensing: strong - galaxies: elliptical and lenticular, cD - 
galaxies: haloes - galaxies: individual: SDSSJ 1430+4105 



1 INTRODUCTION 

Early-type galaxies contain a large fraction of the total 
stellar mass observed in the Universe (e.g., Fukugita et al. 
1998; Bell et al. 2003). Studying the internal structure of 
early-type galaxies is crucial for understanding the baryonic 
physics that plays a key role in the formation and evolution 
of these objects. Several studies have shown that the stars 



assembled in early-type galaxies are embedded in massive 
dark matter haloes (e.g.,Gavazzi et al. 2008; Lagatutta et 
al. 2010; Weijmans et al. 2008), but the precise amount of 
dark matter contained in the galaxies' inner regions is still 
under debate. 

Dark matter only simulations have found indications of a 
universal density profile for dark matter haloes, present 
also in galaxies (the so called NFW profile; Navarro et al. 
1996). Nevertheless, more recent and realistic simulations 
that include also the physics of baryons (e.g., Gnedin et al. 
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2004; Duffy et al. 2010; El-Zant et al. 2001; Bertin et al. 
2003; Ma & Boylan-Kolchin 2004; Blumenthal et al. 1986; 
Jesseit et al. 2002), like radiative cooling and supernova 
and black hole feedback, have suggested that the inner 
profile of the dark matter component can be significantly 
affected by the interactions between baryonic and dark 
matter. 

The internal structure of nearby early-type galaxies has 
been for decades the object of intense dynamical analyses 
(e.g., Saglia et al. 1992; Gerhard et al. 2001; Thomas et al. 
2007; Thomas et al. 2009; Thomas et al. 2011; Pu et al. 
2010). One focus of these studies is to compare stellar with 
dynamical mass to light ratios. The dynamical studies, e.g. 
Gerhard et al. (2001), Thomas et al. (2011) find ratios 
for nearby elliptical galaxies of M/Lb ~ 4 to 10. Similar 
values are also found by Cappellari et al. (2006). Only in 
the last few years has gravitational lensing also contributed 
significantly to our understanding of the luminous and dark 
matter composition of early-type galaxies beyond the local 
Universe (Barnabe et al. 2009; Grillo et al. 2010; Barnabe 
et al. 2010). Strong gravitational lensing in early-type 
galaxies has also proved to be a powerful cosmological tool 
to probe the geometry of the universe independently from 
other diagnostics (e.g., Grillo et al. 2008a; Suyu et al. 
2009, 2010). 

By combining strong gravitational lensing and stellar 
dynamics in a sample of first 15, then 58 early-type 
galaxies, Koopmans et al. (2006) and Koopmans et al. 
(2009) have found that the average total (luminous and 
dark) mass density distribution within the effective radius 
- the radius of the isophote containing half of the total 
light of the galaxy - is well represented by an isothermal 
distribution (p oc r~ 2 ), although significant deviations 
from this result can be observed in individual galaxies. 
Only rare systems where an extended or several distinct 
sources are gravitationally lensed over an extended radial 
area on the lens plane can be used to determine the total 
mass density profile of the lens galaxy over a wide radial 
range through lensing only (e.g. Grillo et al. 2008b, 2010; 
Fadely et al. 2010). Moreover, combining lensing total mass 
measurements with photometric stellar mass estimates 
in these systems offers a unique way to disentangle their 
luminous and dark matter components. 
In this paper, we study the gravitational lensing system 
SDSS J1430+4105 that is composed of a massive early-type 
galaxy acting as a lens for an irregular background source. 
This galaxy was part of the SLACS survey 1 and has been 
studied as part of their lens sample, especially in Bolton et 
al. (2008), Auger et al. (2009) and Auger et al. (2010): 
Bolton et al. (2008) fit a singular isothermal lens model 
to the observed multiple images, while Auger et al. (2009) 
and Auger et al. (2010) combine the measured Einstein 
radii and masses with photometric and dynamical data. 
The surface brightness distribution of the lensed source 
shows several peaks that extend from 4 to 10 kpc from 
the lens galaxy centre. This fact provides the opportunity 
to investigate the lens galaxy mass distribution on radial 
ranges larger than those explored in similar analyses of 
other gravitational lensing systems (e.g., Xanthopoulos et 
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al. 1998, Cohn et al. 2001, Grillo et al. 2010). 
The paper is organised as follows: Section 2 gives an 
overview of the observations and data used in this work 
and introduces the environment of SDSSJ 1430+4105; 
Section 3 describes the details of the strong lensing models. 
Section 4 states the implications for the total mass and the 
mass profile. In Section 5 the mass to light ratio for the 
de Vaucouleurs component is constrained, in Section 6 the 
results are discussed. Appendix A contains further variants 
of strong lensing models. Appendix B gives further details 
of the environment implementation. 

The cosmological model adopted in this paper 
is parametrized by Qm = 0.3, Qa = 0.7, Ho = 
70kms~ 1 Mpc -1 . In the cosmology assumed, 1 arcsec 
in the lens (zi = 0.285) and source (z s = 0.575) plane 
corresponds to 4.30 and 6.55 kpc. 



2 OBSERVATIONS 

The SLACS Survey 2 aimed at finding strong gravitational 
lenses among the galaxies observed in the SDSS. The lens 
detection strategy is presented in Bolton et al. (2004) and 
is based on the examination of the SDSS galaxy spectra, 
taken with a 3 "diameter fibre, to identify emission lines 
not associated with the primary target galaxy but with an 
additional source, aligned with the first galaxy and located 
at a higher redshift. The lens candidates are then ranked in 
terms of their probability of being lensing systems and are 
consequently observed with the HST/ACS and WFPC2. 

Up to now, 85 confirmed (grade-A) lenses (Bolton et 
al. 2006, Auger et al. 2009) were discovered in this way, 
and SDSSJ 1430+4105 is one of these. In Fig. 1 we show the 
SDSS spectrum, from which lens and source redshifts of z\ = 
0.285 and z s = 0.575 are measured, together with the lens 
aperture velocity dispersion of usdss = (322 ± 32) kms -1 . 

2.1 Galaxy light profile and lensing observables 

The basic photometric and spectroscopic properties of 
SDSSJ 1430+4105, taken from Bolton et al. (2008), arc 
stated in Table 1. For these, Bolton et al. (2008) fitted a 
de Vaucouleurs (de Vaucouleurs 1948) profile with ellip- 
tical isophotes to the galaxy's surface brightness distribu- 
tion. They obtained an effective radius of 9 e s = 2.55" = 
10.96 kpc, a minor to major axis ratio of = 0.79, and a 
major axis angle of O q ,L = —12.8°. The angles are trans- 
formed to the adopted local reference frame shown in Fig. 
2, and measured counterclockwise with the y-axis equals to 
0°. 

We retrieve the public HST images from the Hubble Space 
Telescope archive at ESO . Three filters were available 
for this system: HST/WFPC2 F606W with a total inte- 
gration time of 4400s (uall4501m, uall4502m, uall4503m, 
uall4504m), HST/ACS F814W with a total integration time 
of 2128s (j9op36010) and HST/WFC3 F160W with a total 
integration time of 2497s. For the lensing analysis we use 

2 http://www.slacs.org 

3 http://archive.eso.org/archive/hst/ 
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Figure 1. The SDSSJ 1430+4105 spectrum as observed by the 
SDSS 3 "diameter fibre. On the bottom, the observed wavelength 
is stated, while on the top, this is converted into the restframe 
wavelength of the lens (z = 0.285). The dotted vertical lines give 
the SDSS emission / absorption line sample at the redshift of the 
lens. The dashed vertical lines give some prominent emission lines 
at the redshift of the source. Overplotted in the lower part of the 
figure is the flux uncertainty given again by SDSS. The spectrum 
shows several absorption lines typical of an early type galaxy at 
z = 0.285 and some additional emission lines at redshift z = 0.575 
(e.g. the lines at 5872A, 766lA and 7813A, which can be identified 
as the redshifted [OII]3728, H/3 and [OIII]4960 lines respectively). 
Data taken from the Sloan digital sky survey, www.sdss.org. 



the ACS F814W filter observations, since the PSF of the 
ACS camera is smaller than the one of the WFPC2 and 
WFC3. First, we subtract the lensing galaxy's light con- 
tribution with GALFIT (Peng et al. 2002) by using a de 
Vaucouleurs profile, with the parameters of Table 1. Then, 
in order to refine the lens galaxy subtraction and especially 
remove the residuals still present in the central region, an 
additional Sersic profile (Sersic 1963) with index 1.2 is sub- 
sequently subtracted. 

Fig. 2 shows the final galaxy subtracted image. The lensed 
source has a complex surface brightness distribution, with 5 
surface brightness maxima which are imaged 2 times each. 
We mark and label the 5x2 multiple image positions, iden- 
tified as the brightest pixels, in Fig. 2. Their coordinates are 
reported in Table 2 together with approximate error esti- 
mates. We assume in the following that all subcomponents 
A-E are at the same redshift and not unlikely line-of-sight 
projections at different redshifts. The distances of the mul- 
tiple images from the centre of the lens galaxy light dis- 
tribution span a range from 0.93" to 2.32". In the rest of 
the paper, if not otherwise stated, we adopt the coordinate 
system introduced in Fig. 2 which is rotated relative to the 
WCS J2000 (world coordinate system) by 47.21°. 



2.2 Observed environment 

SDSSJ 1430+4105 is not an isolated galaxy. It coincides in 
redshift and location with a galaxy group at z = 0.287, listed 
in the maxBCG cluster catalogue, (Koester et al. 2007). 
Therefore, we should consider the light deflection by the 
lens' environment when we model this lens. We show the 
environment of SDSSJ 1430+4105 (labelled as A) in Fig. 
3. The galaxy labelled as I (and called GI in the following) 
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Figure 2. The multiple image systems that are identified after 
the lens galaxy subtraction and used as input for the lensing anal- 
ysis. The labels Al to E2 mark the positions used for the strong 
lensing analysis. The same letters correspond to images coming 
from the same source feature, the labelling is done according to 
Table 2. The cross marks the centre of the subtracted lens galaxy 
light. Also indicated are the derived shear direction from Sec. Bl 
and the direction of galaxy I. For orientation north is given as 
well. Angles are measured in the local coordinate system coun- 
terclockwise as (-x) over y if not otherwise stated. The image is 
rotated relative to the WCS J2000 by 47.21°. 



Table 2. Observed positions of the multiple image systems 



ID 


01 1 


e-2 1 


Zs 


5 e 


d 1 




(") 


(") 




(") 


(") 


Al 


-1.99 


-0.32 


0.575 


0.05 


2.02 


A2 


0.69 


0.62 


0.575 


0.05 


0.93 


Bl 


-2.08 


0.47 


0.575 


0.05 


2.13 


B2 


1.08 


0.08 


0.575 


0.05 


1.08 


CI 


-2.28 


0.42 


0.575 


0.05 


2.32 


C2 


0.93 


0.03 


0.575 


0.05 


0.93 


Dl 


-1.84 


0.27 


0.575 


0.05 


1.86 


D2 


0.84 


0.80 


0.575 


0.05 


1.35 


El 


0.39 


-1.21 


0.575 


0.05 


1.27 


E2 


-1.64 


1.11 


0.575 


0.05 


1.98 



1 relative to the centre of the galaxy light distribution 

was proposed to be the brightest cluster galaxy (BCG) of 
this group found in Koester et al. (2007). The photometric 
redshift of the group GI is estimated to be z — 0.287 with 
a typical redshift error in the maxBCG catalogue of 0.01. 
Within this error the photometric redshift of the group is 
identical to the spectroscopic redshift (0.28496) of the lens- 
ing galaxy. The group consists of 12 members within the 
R200 of this group, N( gal 2 oo) = 12. At this richness level, the 
maxBCG cluster is typically more than 90% pure and com- 
plete, based on tests with mock catalogues. We now estimate 
the group members based on astrometric and photometric 
data from SDSS DR7 (Abazajian et al. 2009). We consider 
each galaxy within 10' from the main lens J1430+4105. We 
allow for galaxies which have at least one of the photometric 
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Table 1. Photometric and spectroscopic quantities of the lens system 



RA 


Dec 






»?eff 


°"SDSS 


(J2000) 


(J2000) 


(") 


(°) 


(") 


(km s — 1 ) 


14:30:04.10 


+41:05:57.1 


0.285 0.575 0.79 


-12.8 1 


2.55 


322 ± 32 



Given are the position of the galaxy(RA, Dec), the redshifts of galaxy and 
source (z\ z s ), the axis ratio(gL)> the orientation (0 q ,L)i the effective radius 
($ e fT ) of the lens' light distribution and the velocity dispersion <rgrjss- 
Values are taken from Bolton et al. (2008) 
lr This angle is equivalent to —59.3° in the WCS coordinate system, defined 

as (-E) over N. 



redshift estimates (template based (Template-z) (Adelman- 
McCarthy et al. 2007) or neural network based (CC2 z 
and Dl z) (Oyaizu et al. 2008)) consistent within one stan- 
dard deviation with the spectroscopic redshift value of SDSS 
J1430+4105 and the photometric redshift of the group GI. 
The neighbours which pass these requirements are listed in 
Table 3. This table shows that for our definition the galaxy 
A is a group member from both its photometric and spec- 
troscopic redshift. The magnitudes of the galaxies A and I 
in the ugriz filters are 20.43, 19.02, 17.74, 17.12, 16.87 and 
22.54, 19.44, 17.92, 17.35, 17.00 for A and I respectively, 
and thus A is formally the brightest galaxy ('BCG') of this 
group. Since the group membership and group redshift es- 
timate of the maxBCG catalogue (Koester et al. 2007) is 
mostly based on the g-r colour, the contaminated g-r colour 
of A due to lensing has likely led to A not being consid- 
ered as a group member and therefore as the BCG. Using 
the same definition for the r2oo as stated in Koester et al. 
(2007) we find 11 group members from Table 3. 
From Johnston et al. (2007), who tested the maxBCG clus- 
ter finder on simulated groups and clusters of galaxies, we 
derive the probability that a group of the richness given 
in the maxBCG catalogue N( gal:2 oo) — 12 is centred on 
the correct BCG in the maxBCG cluster catalogue to be 
p c (N 2 oo = 12) = 0.63. Therefore both A or I could be the 
true mass centre of the group GI. Finally, the Koester et 
al. (2007) group catalogues could also contain false pos- 
itive detections. Song et al. (2012) have shown that the 
false detection rate of such groups can be as large as 40%. 
This motivates why we will consider lens models with and 
without an external group contribution. We will show that 
including a galaxy group centred on galaxy I has only minor 
influence on the lens parameters. 



3 STRONG GRAVITATIONAL LENSING 

In this section we model the lens mass distribution with the 
public GRAVLENS (Keeton (2001a)) code (Sec. 3.1) assum- 
ing point sources and with the lensview (Wayth & Web- 
ster (2006)) code (Sec. 3.3) using the 2-dimensional surface 
brightness distribution of the same system. Both approaches 
give consistent results. We give a description of the influence 
of the environment on the lens model of SDSSJ 1430+4105 
in Sec. 3.2. 



Figure 3. The observed environment of J1430+4105, taken from 
SDSS. All galaxies with a photometric redshift consistent with the 
spectroscopic redshift of SDSSJ 1430+4105 (labelled with A and 
marked with a diamond) within 1 standard deviation are marked 
by squares. The galaxy marked as I is the BCG for the group 
identified by Koester et al. (2007). In the image, north is up and 
east is left. This image is rotated relative to Fig. 2 by 47.21°. The 
size is 7.5' X 7.5'. For the properties of the galaxies, see Table 3. 



3.1 Parametric modelling using GRAVLENS 

GRAVLENS 4 is a public code that uses parametric lens mod- 
els to reconstruct the properties of an observed lensing sys- 
tem. The lens-modelling we implement here is similar to 
the one of Grillo et al. (2010) where the reader can find 
more details. In this subsection, we use peaks in the sur- 
face brightness distribution of the lensed images as point- 
like position constraints for the lens model (see Table 2 and 
Fig. 2). Since the complex surface brightness distribution of 
the lensed galaxy makes it difficult to associate reliably a 
flux measurement to each multiple image, we neglect flux 
constraints. In GRAVLENS the convergence k for a singular 
isothermal ellipsoid (SIE) or an ellipsoidal powerlaw (PL) is 
parametrized as 



see http:/ /redfive. rutgers.edu/~keeton/gravlens 
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Table 3. Observed environment of SDSSJ 1430+4105 



Obj 


RA 

(J2000) 


Dec 
(J2000) 


magi 


<7i 


Template z 


Error 


NN CC2 z 


Error 


NN Dl z 


Error 


A 


217.51704 


41 


0992 


17.216 


0.006 





228 


0.033 


0.231 


0.027 


0.234 


0.026 


B 


217.44551 


41 


0456 


18.729 


0.018 





202 


0.015 


0.216 


0.073 


0.260 


0.067 


C 


217.45406 


41 


0460 


18.495 


0.020 





262 


0.019 


0.287 


0.045 


0.316 


0.048 


D 


217.45941 


41 


0625 


18.657 


0.019 





179 


0.032 


0.312 


0.055 


0.309 


0.064 


E 


217.46178 


41 


1010 


18.870 


0.016 





266 


0.019 


0.287 


0.055 


0.343 


0.073 


F 


217.46445 


41 


0738 


19.134 


0.020 





286 


0.082 


0.309 


0.096 


0.344 


0.073 


G 


217.47010 


41 


1054 


18.708 


0.013 





224 


0.011 


0.228 


0.073 


0.267 


0.063 


H 


217.49009 


41 


0825 


19.356 


0.024 





200 


0.024 


0.259 


0.070 


0.339 


0.060 


I 


217.49493 


41 


1044 


17.348 


0.008 





273 


0.019 


0.272 


0.020 


0.272 


0.016 


.1 


217.50190 


41 


0995 


18.697 


0.018 





251 


0.021 


0.293 


0.062 


0.326 


0.063 


K 


217.50545 


41 


0779 


18.283 


0.015 





263 


0.020 


0.267 


0.037 


0.293 


0.032 


L 


217.51212 


41 


0655 


18.008 


0.0100 





268 


0.015 


0.273 


0.022 


0.279 


0.023 


M 


217.52402 


41 


1253 


19.358 


0.026 





298 


0.084 


0.451 


0.042 


0.424 


0.042 


N 


217.54462 


41 


1407 


18.947 


0.026 





275 


0.063 


0.284 


0.057 


0.296 


0.048 


O 


217.54641 


41 


0450 


18.500 


0.016 





169 


0.056 


0.210 


0.043 


0.236 


0.057 


P 


217.58642 


41 


0746 


19.373 


0.027 





371 


0.088 


0.401 


0.059 


0.403 


0.064 



Spectro-z Error 



0.28496 



2.13E-4 



Properties of the galaxies considered as part of the environment of J1430+4105. Given in the rows are the object name (Obj), the position (RA, 
Dec), its i band magnitude mag; and its error a\, and its various photometric redshift estimates together with its errors: first the template method 
and its error, summarised in Adclman-McCarthy et al. (2007), followed by the neural network based methods CC2 and its error and Dl and its 
error (Oyaizu et al. 2008). 



K(6i,e 2 ) = 



2(i-)^(S + ei + ffy 



rr (1) 



with 



1 + q 2 ' 

where b is the lensing strength, /3 denotes the steepness 
of the density profile (/? = 2 in the case of an isothermal 
profile), Oi and 02 are the coordinates on the plane of the 
sky relative to the centre of mass of the lens, 9 C is the core 
radius and q is the axis ratio of the isocontours of the con- 
vergence (q = 1 for a circular mass model). In the special 
case of a circular lens without core radius, b equals the Ein- 
stein radius OEm of the lens defined as k(0 < 6>Em) = 1. 
Further we use a de Vaucouleurs Model (de Vaucouleurs 
1948) parametrized as 

I(R) = Iee- 767 [ ( 4F) 1/4 - 1 ] , (2) 

with R cf f being the effective radius (the radius which con- 
tains half the light) and I c the surface density at this radius. 
In GRAVLENS this is implemented as 

1/4 



bdeve 



2)1/! 



©eff 



(3) 



In this parametrisation, bd c v is the value of the central 
convergence. The Einstein radius, however, depends also on 
6 e fi and q. Also a Navarro, Frenk and White (NFW) profile 
(Navarro et al. 1996) is used, defined as 



P(r) 



<5cpc 



r/r a (l + r/r s ) 2 



(4) 



with p c denoting the critical density of the universe at 
the redshift of the lens, and r s and S c are characteristic 
properties of the individual halo. For an overview of its 



lensing properties, see Wright & Brainerd (2000). 

The relation for the LOS projected surface mass density 
£ of the lens and lensing convergence k is 



with 



^crit — 



4ttG Dd-Dds 
c 2 D s 



where Dd, D s and Dd s are the angular diameter distances 
from the observer to the lens, the source and from the lens to 
the source, respectively. The goodness of a model is judged 
by the \ 2 - 



2 

Xlcns 



= E 



Qi-e,u || 2 



(5) 



where ©i denote the model-predicted positions of the i- 
th images, ©04 is its observed position, and &e i its observed 
positional uncertainty. 

Any priors described in the text are added to this x 2 m 
Gravlens via 



with 



2 _ 2 1 2 

Xtot Xlens T Xpr 



2 _ (P ~ Pprior) 

Xprior 2 



where p is the used parameter value, p pr i or its prior 
and ffprior its la error. Results give best-fitting parameters 

and their la errors. The likelihood of a parameter set is 

— 2 11 

given by L oc e Xtot/ . In almost all cases, the best-fitting 
values are within the 68% error interval of the marginalised 
distributions. 

The values of the parameters for the minimum x 2 mod- 
els are given in Tables 4 and 5. There, we give the model 
number, type, the best-fitting parameters of the model to- 
gether with the resultant \ 2 > the number of degrees of free- 

9 2 

dom (d.o.f.) and the reduced x rc d = d o f °^ eacn model. 
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Also, the la error intervals are given. These error esti- 
mates of the parameters are carried out using Monte Carlo 
Markov Chains (MCMC) with several thousand steps each. 
For each model, 10 chains arc calculated with different start- 
ing points. Convergence is reached by comparing the vari- 
ance of the point distribution of each of this chains with its 
combined distribution, see Fadely et al. (2010), Gelman et 
al. (1995). From the final chains, the 2nd half of each chain 
is combined to the final MCMC point distribution. The ac- 
ceptance rate typically lies between 0.25 and 0.3. We explore 
potential parameter correlations from these and derive 68% 
confidence intervals 5 on the parameters by exclusion of the 
lowest and highest 16 % of the MCMC points' distribution; 
the central value is given by the median value of the MCMC 
points' distribution, since there are only small deviations be- 
tween the median and the average values of the 68% and 90% 
error intervals. 

We describe the most important different models without 
environment in the following: To check for the basic proper- 
ties of the system, we model the lens as one component SIE 
(Model I) and PL (Model II) model. To derive the de Vau- 
couleurs masses in this lens, we combine a de Vaucouleurs 
component with a dark matter halo model (Model III) and 
show that this result is not significantly affected by also tak- 
ing the environment into account (Models IV, V). 

Model I : The lens is modelled as a SIE (Eq. 1 with fS = 2); 
the environment of the lens is ignored. The free parameters 
of this model are the lensing strength b, the axis ratio q, 
and its position angle 6 q . The best-fitting model is shown 
in Fig. 4. The results of the MCMC are shown in Fig. 5. 
The density contours describe the probability density for the 
parameter values, whereas the best-fitting model is marked 
with a cross. The reason for the apparent correlation of q 
and b in Fig. 5 lies in the definition of k in Eq. 1. The 
marginalised 68% confidence errors are: b = (1.49^q'q2)"> 
q = (0.71±2;21) and 9 q = {-21.8^2%)° ■ These values are 
in very good (« la) agreement with the values derived by 
Bolton et al. (2008) using a similar parameterisation for the 
lens total mass distribution. 

Model II : Model II follows a power law (PL) (Eq. 1 with 
arbitrary /3 within the limits [1, 2.7]), and thus has one more 
free parameter relative to Model I. The values for the pa- 
rameter distributions are shown in Fig. 6. The marginalised 
distributions change to b = (2.12±J5'.m)", q = (0.81± ?), 
9 q = (-22.21^)° and /3 = (1.73t° il). We observe again 
(see Fig. 6) that the parameters 6, q and f3 are correlated 
with each other. This is entailed by the definition of the 
convergence k in Eq. 1. The steepness parameter /3 is con- 
strained to a value shallower than isothermal on a 1.3a level. 
The orientation O q stays at the same angle as in the SIE 
case, while its axis ratio moves towards rounder solutions, 
now being comparable to the axis ratio of the light distribu- 
tion. 

Model III : In the following, we split the mass distribu- 
tion into different components. We use a de Vaucouleurs 

5 all given errors in this section are the 68% confidence values of 
the marginalised distributions, unless otherwise stated 




Figure 4. The best-fitting model derived with GRAVLENS for the 
SIE case (Model I) (see Table 4), both the image plane and the 
source plane are plotted. On the image plane triangles mark the 
input positions, while crosses mark the best-fitting model posi- 
tions. On the source plane, the predicted source positions are 
plotted. Further, the respective critical lines (caustics) on the im- 
age (source) plane are plotted. 



like component as traced by the stellar component and add 
dark matter with different profiles if needed. Since the de 
Vaucouleurs component for galaxy A alone does not pro- 
vide a good model, see Appendix A, we add a dark mat- 
ter component centred at galaxy A. We add an elliptical 
NFW-likc component to the de Vaucouleurs profile. This 
composition resembles the common picture of galaxy mass 
distribution. For the dark matter halo, we impose a prior 
on the axis ratio based on the Bolton et al. (2008) work of 
qdark, prior = (0.79 ± 0.12). Also we set the limit of the scale 
radius to values < 500", approximately 10 times the value 
we find from Sec B2 for the scale radius. The total mass of 
the de Vaucouleurs component is Md c v = (8.8^'g) x lO n M0 
while the parameters of the dark matter halo are, see Fig. 
7: q d = (0.72±g;i), 9 q = (-26.0±2:i)°, c d = (1.8tJ;°) and 
1200 = (406t 1 29)"- We note that there 

is some degeneracy 

between the concentration c and r2oo- Further we have no 
constraints on r2oo from the data, since we do not have ob- 
servables outside 2.32". Using a NSIE-like dark matter com- 
ponent yields similar results, as described in Appendix A, 
Model Illb. 



3.2 Lens modelling of the environment 

As mentioned before, this galaxy is not an isolated field 
galaxy, hence we investigate the possible impact on the de- 
rived lens parameters by taking the environment into ac- 
count. In the following, we centre a smooth group contribu- 
tion at galaxy I and calculate its convergence and shear at 
the position of SDSSJ 1430+4105. Further modelling of the 
group contribution by summing up the contributions of the 
individual members ("clumpy group"), and by centreing it 
at galaxy A itself is discussed in Appendix B. 

Smooth group mass distribution centred at galaxy I 

According to Rozo et al. (2009), we can transform the group 
richness into a group mass of M 5 oo = (0.72 ±0.29) x 1O 14 M 
within la. This mass can be converted into a velocity disper- 
sion of crgroup = 519 ± 107kms~ 1 , using the critical density 
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Table 4. minimum-x 2 values and parameter estimates derived with GRAVLENS for the isothermal and powerlaw models 
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Table 5. minimum-x 2 values and parameter estimates derived with GRAVLENS for the two component de Vaucouleurs plus dark matter models 
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Figure 5. Error estimates of the MCMC for the SIE case (Model 
I), plotted are the individual points of the MCMC together with 
the 68 and 90 % confidence regions of the distribution, indicated 
by the density contours. The crosses mark the minimum-x 2 value 
from Table 4. The bars on the axes mark the respective 68 % 
marginalised error intervals. 



of the universe 

con i \ 3M 50 o 

47rr 500 

and the singular isothermal sphere equation: 




0.8 



0.7 




-30 -25 -20 

©,0 





1.6 1.8 2.0 2.2 



1.6 1.8 2.0 2.2 



Figure 6. Error estimates of the MCMC for the PL case (Model 
II), plotted are the individual points of the MCMC together with 
the 68 and 90 % confidence regions of the distribution. The crosses 
mark the minimum-x 2 value from Table 4. The bars on the axes 
mark the respective 68 % marginalised error intervals. 
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Figure 7. Error estimates of the MCMC for the NFW+de Vaucouleurs model (Model III), plotted are the 68 and 90 % confidence regions 
of the distribution. The crosses mark the minimum-x 2 value from Table 5. The bars on the axes mark the respective 68 % marginalised 
error intervals. The individual points of the MCMC arc omitted for clarity. 



There p c (z) denotes the critical density of the universe at 
redshift z and a gTOUp the velocity dispersion of the group. 
Subsequently this gives an Einstein radius of Oe = 3.6 ± 
1.5", using again a SIS assumption (see Sec. 3.1 for details). 
This results in a convergence of and shear of 



SIS group 
A 



0.029 
7 sis group = 002 9 



(6) 



al. (2007). We obtain a convergence and shear of 



^NFW group 



NFW group 

7a 



= 0.025 
= 0.026 



(7) 



Further, we note that the angle of A towards I is —26°, 
therefore forming an angle of 16° with the external shear 
value derived in Sec. Bl. We examine the HST and SDSS 
frames which cover galaxy I and its vincinity for group- 
scale multiple images to further constrain the group mass 
distribution but do not find any sign for strong lensing. 



at galaxy A if galaxy I is assumed to be the group centre. 
Alternatively, we model the smooth group as a 'typical' rich- 
ness 12 galaxy group NFW (Navarro et al. (1996)) halo with 
concentration c = 4.22 and r2oo = 848 kpc from Johnston et 



Model IV : From Section 2.2 we expect that there is some 
environment dark matter present in this galaxy. We check 
whether using this group dark matter contribution with a 
de Vaucouleurs component for galaxy A is sufficient to ex- 
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Table 6. minimum-x 2 values derived with LENSVIEW 




6 7 8 9 

bgroup ( ) 

Figure 8. Error estimates of the MCMC for the de Vau- 
coulcurs+GI model (Model IV), plotted are the individual points 
of the MCMC together with the 68 and 90 % confidence regions 
of the distribution, indicated by the density contours. The crosses 
mark the minimum-x 2 value from Table 5. The bars on the axes 
mark the respective 68 % marginalised error intervals. 

plain the observations, even though modelling this system 
with a pure de Vaucouleurs component fails, see Appendix 
A. Therefore, in this model, we combine the de Vaucouleurs 
profile with a group halo centred at galaxy I. To account for 
the environment, we include the galaxy group explicitly as 
a SIS profile centred at galaxy I in Table 3. We use a prior 
on the group Einstein radius of b gI . ouPi p rio r = (3.6 ± 1.5)". 
The de Vaucouleurs component has shape parameters as 
stated in Table 1. The group acts almost as a mass sheet. 
We get a x 2 — 18-9 for the best-fitting model. We get 
parameter estimates of M dc v = (13-51Q?,) x 1O 11 M0 and 
bgroup = (8.OI0 7)" as can be seen in Fig. 8. Besides be- 
ing a worse fit than most of the other models, this model 
also needs a much more massive group present than what is 
likely from the observations. Therefore, dark matter that is 
distributed almost uniformly within Ge of the galaxy does 
not provide a good model for the system. 

Model V : This model adds environmental effects to Model 
III. Therefore we add the group GI explicitly as above, yield- 
ing 3 components: the group GI, the dark matter associated 
with the galaxy as an elliptical NFW profile and a stellar 
component modelled as a de Vaucouleurs profile. We use 
the same constraints as for Model III. We get the following 
parameters, see also Fig. 9: M dc v = (10. 4^^) x 10 11 M o , 
q d = (0.7912:11), 9 q = (-21.61?;!)°, c d = r 200 = 

(32111^)", and for the galaxy group b group = (4.911;?)". 
We note that these parameter estimates do not significantly 
change compared to Model III, therefore the inclusion of 
group GI has only a small influence on the estimated galaxy 
parameters; the M de v for the de Vaucouleurs component 
is slightly increased. Again, we are not able to constrain 
the concentration c or r2oo of the dark matter halo. Mod- 
els Va and Vb in Appendix A employ a NSIE-like galaxy 
dark matter halo (Model Va) and an external shear contri- 
bution instead of a explicit group contribution (Model Vb) 
and again give results very similar to Model V regarding the 
parameters for the lensing galaxy. 

3.3 Full surface brightness distribution using 

LENSVIEW 

We also use lensview (Wayth & Webster (2006)) to derive 
models and mass estimates for SDSSJ 1430+4105 and to 
reproduce the full surface distribution of the lensed galaxy 
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and its unlensed source. Lensview is a publicly available 
program that fits parametric lens models to image data and 
uses the best-fitting lens model to reconstruct the source and 
image. The code uses the image data, a corresponding noise 
map, and an image mask to minimise x 2 — \S, where x 2 is the 
chi square difference between the reconstructed image and 
the data, S is the entropy in the source, and A is internally 
adjusted such that \ 2 approaches its target value. If the data 
are well fit by the model, the entropy term serves to smooth 
the source. Because the flux of each unmasked data pixel 
is used in the fit, lensview is well-suited to systems with 
extended flux like SDSS J1430+4105. The profile used here 
is defined, following Barkana (1998), as 

^e^^^ + f)^, (8) 

where b' gives the Einstein radius, q the axis ratio and /3 
again the power law exponent of the profile. We note that the 
normalisation of the profiles is different from Eq. 1, resulting 
in different values for the Einstein radius in both approaches. 

The minimum-^ 2 results are stated in Table 6. The SIE 
best-fitting parameter values derived here agree with those 
found in Sec. 3.1, when directly compared to Model I in Ta- 
ble 4. For the PL model, we see a consistency of the differ- 
ent models from gravlens and lensview within the stated 
errors for q, Bq and j3. Since, as mentioned before, the nor- 
malisation of the convergence profiles is different, the 6 / 
b' values do not compare directly to each other. For the 
models including the environment, the direct comparison of 
the SIE+ES model with Model la shows again a consistency 
within the errors derived in Appendix A for the lens parame- 
ters. However, the external shear angle shows a discrepancy, 
the angle is offset relative to the expected value derived in 
Sec. 2.2. Since the external shear contributes no mass, this 
will not have a significant effect on the mass estimates in 
Sec. 4. The same is true for the PL+ES case, the compar- 
ison with Model Ha gives an agreement within the given 
errors in all parameters besides 6 7 . 

3.4 Tests on the strong lensing assumptions 

For all but the single, pure de Vaucouleurs model, the centre 
of the light and mass distribution do not necessarily have to 
coincide. To address possible systematic effects of the as- 
sumption that the mass model profiles are centred on the 
centre of light, we use Model I with a free centre with priors 
on the centre position of uncertainty a® 1 1 = ue 2 , = 0.2". 
In this case, we find: b = (1.5llg;{g)", q = (0.69^;^) and 
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Figure 9. Error estimates of the MCMC for the de Vaucouleurs+NFW model (Model V), plotted are the 68 and 90 % confidence regions 
of the distribution. The crosses mark the minimum-x 2 value from Table 5. The bars on the axes mark the respective 68 % marginalised 
error intervals. The individual points of the MCMC arc omitted for clarity. 



O q = (— 26. ll 2 ,^) - The lens total mass centre moves to 
0i,i = (-0.02t°;°*)" and 2 ,, = (-0.111° °*)". The lensing 
strength b and the axis ratio q agree within the errors with 
the values derived for Model I, but the uncertainty on q in- 
creases. While the lens position 0i,i is still consistent with 
0, there is a dependency of 02, i and q . This has no effect 
on the mass profiles M(< R), so we conclude that we can 
fix the lens total mass centre to its centre of light without 
inducing systematic effects on the derived masses. 
In the previous sections we restricted the smooth mass 
model of the group GI to be spherically symmetric. We al- 
ter the shape of GI to an elliptical isothermal profile (SIE) 
and study the changes on the best-fitting SIE lens mod- 
els. This model is analogous to Model la, but combines a 
SIE lens with a SIE instead of a SIS galaxy group centred 



on galaxy I. Therefore, we create 2000 random representa- 
tions of GI as a SIE. For each of these representations, the 
axis ratio is randomly chosen between qGi = [0.6, 1] with 
a random orientation. The centre of this SIE group model 
is taken from a gaussian distribution, centred on galaxy 
I, with a width of <j Xi gi = fy,Gi = 1.5". We recalculate 
the best-fitting parameters for each of these SIE+GI(SIE) 
models and evaluate the scatter of the best-fitting parame- 
ters to estimate the systematic uncertainties emerging from 
the assumption of a simple SIS group halo. This results in 
68% c.l. distributions of the best-fitting parameters for this 
modified Model la: bi ons = [1.44,1.46]", qi cns = [0.79,0.82], 
q = [-18.4,-16.2]° and b group = [4.27,4.42]". These in- 
tervals are small compared to the statistical uncertainties for 
Model la derived from the MCMC sampling. We conclude 
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Figure 10. best-fitting model for the SIE case 10(a), its residuals 
10(b) and the corresponding source 10(c). For the images 10(a) 
and 10(b), the pixel size is the same as in Fig. 2; the source has 
a 3 times smaller pixel size. The same flux scaling has been used 
on all images. 



that the details of the group model representation do not 
change the results for the lensing galaxy significantly, there- 
fore including the most simple SIS model for GI is sufficient. 
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Figure 11. The masses and radial derivatives for the Models 
stated in Table 7 and Table 8. Plotted are the mean values to- 
gether with its 1 a error bars. For Models I and II, the black 
asterixes mark the values at the Einstein radius (fixed to 1.51"), 
while the blue triangles mark the extrapolations to the effective 
radius of the galaxy. For the de Vaucouleurs like models II and 
V, the green squares are the values at the Einstein radius, the red 
crosses are the ones at the effective radius. While the total mass 
within the Einstein radius is tightly constrained independent of 
the model used, the mass within the effective radius depends on 
the mass model used. 



4 RESULTS FOR THE GALAXY MASS 
PROFILE 

The total masses M(< R) within a cylinder of radius R and 
its derivatives obtained from the SL analysis for Models I to 
V are shown in Figs. 12 and 13. We have calculated these 
values within several concentric apertures with radii ([0.46 
0.92 1.26 1.59 1.95 2.32 2.78] "), chosen to lie in the radial 
regions covered by the lensed images plus extrapolations to- 
wards smaller / larger radii. For the Models I to V of Sees. 
3.1 and 3.2, the masses are estimated by randomly taking 
1000 MCMC points and creating convergence maps for each 
one of these 1000 models. The 68 % (90 %) errors are esti- 
mated by taking the central 680 (900) models at each radius. 



4.1 Mass profiles for the single component 
isothermal and powerlaw models 

First we focus on the masses derived for Models I and 
II in Table 4. The Einstein radii are defined as the radii 
within which the mean convergence equals 1. For this, we 
calculate the mean convergence around the Einstein ra- 
dius in 0.03" distance bins. The results of this calculation 
are stated in Table 7. Since all models agree on an Ein- 
stein radius of 0e = 1.51"== 6.48 kpc == Re, we adopt this 



value as "the" Einstein radius of this lens with an uncer- 
tainty of 0.03"= 0.13kpc. We get a mean Einstein mass of 
(5.37 ± 0.06) x 1O U M for the Models I and II with a fixed 
Einstein radius of Oe = 1.51". This values are in good agree- 
ment with the ones stated by Auger et al. (2009) for this 
system also based on strong lensing. 

We also extrapolate the models to the effective radius 
r c fr = 2.55"= 10.96 kpc of the galaxy, and calculate the 
mass and its derivative. We find an enclosed mass between 
Mtot.onc = 8.9 x lO^M© and M to t,enc = 11.3 x 10 n M Q on 
a la level, depending on the model used. The azimuthally 
averaged results of the included masses are plotted in Fig. 
12. For an SIE model, the mass included within radius r 
grows linearly with the radius, so the derivative of it is ex- 
pected to be independent of the radius. This is the case 
for the singular isothermal model (Model I) in Table 4. If 
we allow the steepness to vary (Model II) the mass profile 
tends to be steeper at the Einstein and effective radius. For 
the radial mass derivative at the Einstein radius, we cal- 
culate values between dMt d °^ cnc = 0.8 x lO^Mekpc -1 and 
dMt d °^ cnc = 1.2 x lO^Mokpc" 1 . The extrapolation to the 
effective radius ranges from dMt ^ cnc = 0.8 x lO 11 M kpc~ 1 
for Model I to iMt °^>"= = 1.6 x lO^Mokpc" 1 for Model II. 
This values are plotted in Fig. 11. Here and in Table 4 we 
state the 68 % c.l. errors. 
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Table 7. The masses and mass derivatives at the Einstein radius, the globally adopted Einstein radius and the effective radius for the different 
models. Given are the projected, enclosed masses and its derivatives 
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Figure 12. The mass within radius R, M(< R), from the models I and II, analysed in Sec. 3 is plotted in Fig. 12(a). For each model, 1000 
random entries from the MCMC are used to calculate the errors. Each time, the means and 90 % errors are shown, with the respective 
means in bold lines. The bars at the bottom mark the radii of the apertures used to calculate the enclosed projected masses. The masses 
are in units of 1O 11 M , the radii are stated in kpc in the lens plane. The mass estimates for the Models I and II are shown as solid (I) 
and dashed (II) lines. Vertical lines indicate the Einstein and effective radius in both plots. Fig. 12(b) shows the same as Fig. 12(a) but 
for the LENSVIEW derived masses. The errors are estimated by an increase of the reduced \ 2 of the extended model by 1. As can be seen, 
the masses agree with each other in terms of derived masses within the errors. 



4.2 Mass profiles for the de Vaucouleurs plus 
dark matter halo models 

From the single component lens analyses in Sec. 4.1, we 
conclude that the total projected mass density profile is 
isothermal or slightly shallower than isothermal. The de 
Vaucouleurs mass density drops faster with radius than the 
isothermal profile. Therefore, we expect the pure de Vau- 
couleurs profile to be a poor description of this lens' mass 
profile (as seen in Model IV in Sec. 3.2 and Model IVa in Ap- 
pendix A) and expect that we need to add some component 
which follows a shallower than isothermal density profile. In 
principle a constant mass sheet, like a nearby group halo 
would approximately provide, could do this, given that it 
provides enough mass at the position of the lens. In Model 
IV, we see that this model is a poor fit to the data. Therefore 
a (dark) matter component centred at the position of the 
lens is necessary. We model this in Models III, Ilia and Illb 
in Sec. 3.1 and Appendix A with different types of dark mat- 



ter halo profiles. A SIE-like dark matter component (Model 
Ilia) suppresses the de Vaucouleurs part of the matter pro- 
file, effectively yielding a single component model. Because 
the light distribution is well-fit by a de Vaucouleurs profile, 
we require a non-zero de Vaucouleurs component for this 2 
component fit, hence we do not consider Model Ilia in the 
following. The Models Mb, Va and Vb in Appendix A give 
similar results to Models III and V in Sec. 3.1. Therefore, in 
the following, we mostly consider Models III and V, which 
model the lens using an NFW profile for the dark matter 
component. 

Besides the stars and the dark matter, an elliptical galaxy 
or a galaxy group also contains some amount of gas. Since 
we do not model this component individually, this gas needs 
to be incorporated in either the dark matter or the de Vau- 
couleurs component, effectively limiting the accuracy of our 
mass estimates to the gas mass fraction in elliptical galax- 
ies and groups of galaxies. Young et al. (2011) get typical 
molecular gas masses of early-type galaxies in the ATLAS 3D 
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project of M(H 2 ) < 1O 9 M , less than 1% of the total galaxy 
masses derived here. The hot gas component in an elliptical 
galaxy or a group of galaxies can contribute up to 10% of 
the total mass in the centre of the galaxy or group of galax- 
ies, see e.g. Sanderson et al. (2003). Hence the uncertainty 
of our mass estimates due to the neglected gas is ^ 10%. 
We again adopt 8e = 1.51" for the Einstein radius. First we 
focus on the masses within this radius, see Table 8. For the 
total masses within the Einstein radius of the models III and 
V, we measure M to t,Ein = (5.33 ± 0.04) x 10 n M Q The ra- 
dial mass derivative is 



dM to 



dR 



= 0.86±S;S? x lO^Mokpc- 1 



. The extrapolations to the effective radius give M to t, e ff = 
9.5±°i x lO^Mo for the mass and dMt ° t >— - 1 nn+o.ie v 



1.00: 



dR - —"-0.14 - 

10 M Q kpc for its derivative. This values are plotted in 
Fig. 11. As can be seen, the enclosed masses and its deriva- 
tives at the Einstein radius and the effective radius agree 
with each other throughout Models I, II, III and V. We 
state the de Vaucouleurs mass within the Einstein and ef- 
fective radius of Model III as Component IIIA in Table 8. 
We get a mass of M dc v,Ein = (3.2t°;f) x lO n M , mean- 

M doV ,Ei, 



ing that 



35% of the total de Vaucouleurs mass 



M dcV,tot 

is concentrated within the Einstein radius for this lens. For 
Model V, we get similar values for the de Vaucouleurs com- 
ponent, see Component VA in Table 8. In Fig. 13 and Fig. 
A3 in Appendix A, the projected, enclosed lens masses and 
their derivatives are plotted versus radius for the different 
2 component strong lensing models. The measurements are 
done using circular apertures, so all of these values are az- 
imuthally averaged. As one can see, including an explicit 
group halo at GI (Model V) has only a minor influence on 
the mass estimates and their derivatives. The total masses 
agree very well with the one component estimates in Fig. 12. 
Also, all models agree very well on the total masses and their 
radial derivatives, tending to give a shallower than isother- 
mal mass profile in the centre. For the Models III and V 
in Fig. 13, the dark matter haloes modelled as NFW-haloes 
agree very well with each other, meaning that the environ- 
ment has only minor influence on the mass estimates. This 
is also true for the de Vaucouleurs component. We note that 
the uncertainties on the individual components are larger 
than the uncertainties on the total masses and derivatives, 
giving a well-constrained total mass. 



4.3 3d spherical reconstruction 

Further, we also reconstruct the 3d matter densities from 
the 2d data for Model III. For this, we employ the inverse 
Abel-transform: 



p(r) 



7T Jr 



dT, dR 



dR^W 



(9) 



transforming a 2d circular density function £ into a 3d sper- 
ical density function p. Since this only transforms circular to 
spherical profiles and vice versa, we start from the mass mea- 
surements within a cylinder in Fig 13(a) for the azimuthally 
averaged profile. In Eq. 9 the integration extends to infin- 
ity, which is not possible due to our limited range of reli- 
able data. To estimate the radial range at which we can use 
Eq. 9 only integrating up to our last data bin we test it on 
a SIS toy model. For a SIS toy model, we know both the 
spherical and the projected circular density. We then con- 



sider this radial range reliable where the deviation of the re- 
constructed 3d density from the analytical SIS density does 
not exceed 2 x 10 7 M Q kpc~ 3 . From this comparison, we con- 
clude that this inverse Abel transformation is only reliable 
up to « 6.5 kpc with a systematic error smaller than 30%, 
given our limited radial range of data. The reconstructed 
3d profile is shown in Fig. 14. The errors plotted are only 
statistical, not taking any systematic effects into account. 
The dark matter accounts for only a minor fraction of the 
total mass in the 3d centre of the galaxy. We now turn to 
Fig. 7 in Thomas et al. (2011). In the lower part of Fig. 
7 they have displayed the ratios of the mean dark matter 
density and mean total density within the Einstein radius 
of Coma galaxies as a function of their velocity dispersion. 
(For the definition of the synthetic Einstein radius for Coma 
galaxies, see Thomas et al. 2011). To see whether there are 
structural differences for the Coma and the higher redshift 
SLACS sample one would like to enter the corresponding 
deprojected values for SLACS galaxies in these figures as 
well. These were not available until now because the dark 
to total matter fractions were only calculated for the line of 
sight projected densities within the Einstein radii (i.e. cylin- 
drical averages) by gravitational lensing. The corresponding 
projected values are shown for SLACS (and Coma) galaxies 
in the upper part of Fig. 7 of Thomas et al. (2011). The 
projected and deprojected values differ, since the projection 
along the line of sight mixes scales: Parts of the matter that 
has a large physical distance from the centre of the galaxy 
but lie on the line of sight are taken into account when 
calculating the projected dark matter fractions. Due to the 
monotonic increase of the dark to total matter density ra- 
tios as a function of radius, the projected ratios displayed in 
the upper part of Fig. 7 of Thomas et al. (2011) are upper 
limits to the central, 3-dimensional density ratios at the Ein- 
stein radius. With the analysis described in this work we are 
able to measure the 3-dimensional densities of the (spheri- 
cally averaged) dark matter and de Vaucouleurs components 
of the lensing galaxy separately from gravitational lensing 
alone due to the large radial coverage of multiple images in 
the image plane by one source. Since the source is only one 
background object, we do not need to take the systematic 
uncertainties into account that arise in systems with multi- 
ple image systems from sources at different redshifts like e.g 
Gavazzi et al. (2008). At the Einstein radius we obtain (us- 
ing Fig. 14, displaying Model III (de Vaucouleurs+NFW)) 
a dark to total density ratio of 22 % for the dark to to- 
tal density. Doing the same for the Model V where SDSSJ 
1430+4105 (consisting of de Vaucouleurs and dark matter 
component) is embedded in a DM halo centred on galaxy 
I we find that the ratio of dark to total matter density at 
the Einstein radius is about 14%. Since the dark matter 
fraction increases towards the outskirts, these ratios of den- 
sities at the Einstein radius are upper limits for the mean 
dark matter to total matter density ratios of galaxy SDSSJ 
1430+4105 within the same Einstein radius. On a (90% c.l.) 
basis, the density ratios at the Einstein radius are larger 
than 15% (Model III) and 5 % (Model V) for the dark to 
total matter density. 
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Table 8. Same as Table 7 for the 2 component models. 
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Figure 14. The 3d densities for the stellar (dashed) and the dark 
(dash-dotted) component, as well as the total matter (solid) For 
this, the data of Model III is used. Plotted are the median (thick) 
and the 90 % c.l. intervals. Spherical symmetry is assumed for 
this reconstruction. 



5 MASS TO LIGHT RATIOS FOR THE DE 

VAUCOULEURS COMPONENT AND DARK 
TO TOTAL MASS RATIO 

Since we calculated the de Vaucouleurs masses for this 
galaxy, we now want to estimate the rest frame mass to 
light ratios of this galaxy. Further, we evolve this mass to 
light ratios to present day values in order to compare it with 
Coma galaxies. First, we calculate the dark matter fractions 
within the Einstein radius. In Fig. 15, we plot the dark over 
total enclosed mass fraction within the Einstein radius. The 
error bars are estimated as before from the central 68 % and 
90 % entries of the random sample drawn from the MCMC 
for Models III and V, respectively. The fractions for Model 
III and V are: ^™ = (0.40t^) and ^hew_ = ( .27l^?). 



M t „t - V.""™-0.09^ M tot 

These fractions indicate a substantial amount of dark mat- 
ter within the Einstein radius of this lens. As can be seen, 
this picture is not significantly altered by locating the group 
at GI in Model V. Although the actual numbers change, we 
still need dark matter associated with the lens in the centre 



of the galaxy. Since we ignore the dark matter contribution 
associated with the group GI in Model V for the total mass, 
we get a lower dark matter fraction for Model V compared 
to Model III. 

For the Models III and V, we calculate the mass to light 
ratios for the de Vaucouleurs component at the redshift of 
the lens. We use the masses from lensing and the light (in 
rest frame B and R) as obtained from photometric data. To 
compare the mass to light ratios with present day galax- 
ies, we also need the luminosity evolution to redshift zero in 
these bands. We take the observed griz SDSS photometry 
for this system which covers the rest frame B and R fil- 
ters to calculate the B-R rest frame colour. 6 We calculate a 
colour of (B — R) = (0.80 ± 0.03) and luminosities of L B , r f = 
(1.66 ±0.03) x 10 11 Lq and L R , rf = (1.92 ±0.02) x IO^Lr.q 
from the absolute rest frame magnitudes. To estimate the 
luminosity evolution until today, we fit 3 extinction-free 
BC (Bruzual & Chariot 2003) composite stellar popula- 
tion (CSP) models to the observed griz SDSS photometry: 
A Salpeter initial mass function (IMF) (Salpeter 1955) with 
solar metallicity (Model A) and 2 models with Chabrier IMF 
(Chabrier 2003) and solar / super-solar metallicity (Mod- 
els B,C), respectively, see Drory et al. (2001, 2004). The 
best-fitting results are stated in Table 9 together with the 
luminosity evolutions in B and R bands and the best-fitting 
stellar masses, which agree with Grillo et al. (2009). From 
the best-fitting star formation histories (SFH) to the spec- 
tral energy distribution (SED) we obtain a stellar age of 
typically 8 Gyrs and a B-R colour of 0.8. Therefore, this 
galaxy has a formation redshift of approximately 2 to 3 
which is a typical value for elliptical galaxies. If we divide 
the de Vaucouleurs masses derived by lensing in Sec. 3 by 
the luminosities derived from the SDSS photometry we ob- 
tain the mass to light ratios (MtoL) for the de Vaucouleurs 
component of this galaxy. For Models III and V, we find a 
MtoL of = (5.31°:!)^ and ^ = (6.21°;°)^ 

in the B-band rest frame at the redshift of the lens. These 
two Models give the same MtoL within the errors, although 
including the group GI explicitly increases the most-likcly 
MtoL. We compare this with the total light of the galaxy 
and the stellar mass derived in Grillo et al. (2009), who 
use composite stellar population models with a Salpeter or 
Chabrier IMF, a delayed exponential star formation history, 
and solar metallicity to model the SDSS multi band pho- 



All magnitudes arc given in AB. 
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Figure 13. The mass within radius R, M(< R) from the Models III and V analysed in Sec. 3. For each model, 1000 random entries 
from the MCMC chains are used to calculate the errors. All errors plotted are the 90 % error intervals with the respective means in 
bold symbols. The bars at the bottom mark the radii of the apertures used to calculate the enclosed projected masses or its derivatives. 
The masses are in units of 10 Mq, the radii are stated in kpc in the lens plane. In 13(a), the mass estimates of Model III for the de 
Vaucouleurs (dashed), NFW (dash-dotted) and its sum (solid) are plotted. While the sum of these two is fairly well constrained, the 
errors on the individual parts are bigger. In 13(b), the same mass estimates are plotted for Model V, together with its 90 % error intervals, 
splitted in de Vaucouleurs (dashed), NFW (dash-dotted) and GI (dotted) parts and its total sum (solid). The radial mass derivatives 
are plotted in 13(c) for Model III and 13(d) for Model V, keeping the line coding. Plotted is the change in enclosed mass with radius. 
This can also be interpreted as the mass in a thin ring with width dR at radius R, M(R). Again, vertical lines indicate the Einstein and 
effective radius in both plots. 



tometry. First, we compare in the rest frame B band. In 
Fig. 16 we plot the cumulative distribution function for the 
stellar mass to light ratios derived from the respective de 
Vaucouleurs parts of the Models III and V and Models Illb 
and Va from Appendix A. We overplot the stellar MtoL ra- 
tios derived in Grillo et al. (2009) for this system and get 



the best agreement for a NFW like halo and a Salpeter IMF. 
In the R-band, we get a MtoL for Models III and V of 
= (4.6+°-?) ^ and ^v = (5.4+^)^. To com- 
pare the lensing galaxy SDSSJ 1430+4105 with present day 
Coma galaxies (Thomas et al. 2011) we have to account for 
the luminosity evolution between redshift 0.285 and now. 
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Table 9. Galaxy luminosity evolution factors for the different 
IMF and metallicity models 
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Figure 15. Dark over total enclosed mass ratio for this galaxy. 
Shown are Models III (solid) and V (dashed). For Model V the 
group contribution is ignored, meaning that we only consider the 
dark matter associated with the galaxy itself as dark matter and 
only the sum of de Vaucouleurs and dark matter as the total mat- 
ter. The vertical line indicates the effective radius of the galaxy. 
As one can see, both models predict dark matter contributions 
in the centre of this galaxy, with changing amounts depending 
on the modelling details. Plotted are the 68 % and 90 % errors, 
respectively. 



We use the average evolution factor from the SFH mod- 
els stated in Table 9, derived from the extrapolations of 
the fitted SFH models to redshift 0, which increases the 
MtoL in the R band by a factor of 1.48 for z~0, giving 



/"MdeV \ 



(6.8 



+ 1.2- 
1.6. 



At,: 



L R /prcs-^_i. WL0 ° E <mu l Lr , pr cs - V«-«- 1 .3 ) L R 

for Models III and V, respectively. In Fig. 17, we plot this 
R-band de Vaucouleurs mass to light ratio at redshift zero 
against the present day R-band mass to light ratio for a 
Kroupa-IMF (Kroupa 2001), obtained again from the SFH 
fit of Grillo et al. (2009), translated to R-band and evolved 
to redshift zero. We added the results for a dynamical study 
of Coma galaxies by Thomas et al. (2011). This allows us 
to conclude that SDSSJ 1430+4105 evolves into a galaxy 
with mass to light ratio similar as the Coma galaxies, and 
shows the same conflict with respect to a Kroupa IMF as 
they do. This conflict to a Kroupa IMF would be resolved 
if, for example, the de Vaucouleurs component is not made 
of stars only but contains dark matter as well. 
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Figure 17. MtoL in R from lensing versus stellar mass using a 
Kroupa-IMF, sec Grillo et al. (2009). Both values are evolved 
passively to the redshift 0. The lensing estimate is taken from 
Model III. The black points show the mass to light ratios derived 
from stellar dynamics analysis of the COMA cluster galaxies ver- 
sus their stellar MtoL ratio derived from photometry for a Kroupa 
IMF, taken from Thomas et al. (2011). The solid line marks the 
one to one correspondence. 



6 DISCUSSION AND CONCLUSIONS 

In this paper, we studied the lensing properties of SDSS 
J1430+4105. From the complex source, we identified 5 dou- 
ble image systems, spanning a radial range from below 0.9" 
to almost 2.1". The source is spectroscopically confirmed at 
a redshift of 0.575. Parametric models can match the ob- 
served image positions well with an average scatter in the 
position comparable to the pixel size of the ACS camera 
input image, which is 0.05". Our results are: 

I) The best-fitting reconstruction of the profile favours a 
profile slope shallower than isothermal for the best-fitting 
model. However, profiles with free slope for the density 
steepness are consistent with an isothermal profile at 90% 
c.L. This is also true when combining an explicit model 
for the de Vaucouleurs-like light distribution with a NFW- 
like dark matter component. Auger et al. (2010) found 
a steepness for the 3d density profile for this system of 
p ,sj T .-( 2 -06±o.i8) ky using the location of the Einstein radius 
only and combining this with stellar dynamics, in agreement 
with our results for the one component powerlaw total mass 

-(1 7'j+ - 21 ) 

distribution of p ~ r y -0.131 within the errors. 
II) The galaxy is part of a group of galaxies listed in the 
maxBCG cluster catalogue. Using a lens mass component 
following the stellar light, we can not model the strong lens- 
ing signal for this galaxy if we use this component alone or 
combine it with a dark matter halo not centred on SDSS 
J1430+4105, called galaxy A. Therefore, this leads to two 
possibilities: Either A is indeed the centre of the galaxy 
group, or the galaxy A is a satellite of this group, resid- 
ing in its own dark matter halo. Since we cannot distinguish 
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Figure 16. Cumulative probability distribution functions of the mass to light ratios for the de Vaucouleurs components of the models 
III and V. Models III, Illb, V and Va are marked by red, blue, turquoise and green. The vertical lines mark the derived stellar mass to 
light ratios with its 1 <r errors from Grillo et al. (2009) for this system fitting SFH to broad band SDSS photometry using a Salpeter 
IMF (solid line) and Chabrier IMF (dashed line). The mass to light ratios are as observed at z = 0.285 and not corrected for luminosity 
evolution to redshift zero. 



between these 2 cases, neither from the lensing signal nor 
from external data, we model both scenarios and show that 
these lead to similar results regarding the mass distribution 
of the galaxy. We show that the dark matter halo of galaxy A 
must not be singular and isothermal at the same time, since 
this would suppress the de Vaucouleurs component. Both a 
non-singular, isothermal halo and a NFW-like halo for the 
dark matter halo of galaxy A fit the data well. We find agree- 
ing dark matter fractions and distributions for both cases. 
From the lensing data, we cannot distinguish whether the 
dark matter halo follows a NFW or NSIE profile in the cen- 
tre, since we cannot constrain the concentration c and T2oo 
- or O c for a NSIE dark matter halo - of the dark matter 
component well. From the models taking explicitly the en- 
vironment into account we conclude that the dark matter 
and the light of the galaxy have the same axis ratio and are 
likely coaligned. 

Ill) We estimate the rest frame B-band mass to light ra- 
tios for the lensing galaxy from the de Vaucouleurs lensing 
component. For the case of a deVauc+NFW mass model, 
we obtain a total mass of M dcV ,tot = (8.811;!) x 1O U M 
for the de Vaucouleurs component. Grillo et al. (2009) have 
obtained stellar mass estimates for the luminous component 
using the ugriz broad band SDSS photometry and SFH fits. 
They assumed solar metallicity composite stellar popula- 
tions with a delayed exponential SFH, and examined the 
Salpeter IMF (with BC and MAR (Maraston 2005) single 
stellar populations (SSP)), and the Chabrier and Kroupa 
IMF (based on BC-SSPs). They obtained stellar masses 



of M* = (5.6l?|) x 1O U M , M t = (3.9l^l) x 10 n M Q , 
M* = (3.2t°i) x 1O U M , and M* = (2.91?;®) x IO^Mq for 
these four cases. The stellar mass agrees best with the mass 
of the de Vaucouleurs component obtained from lensing if 
we assume a Salpeter IMF. In principle, the mass to light 
ratios also depend on the age of the galaxy and its metal- 
licity. According to De Lucia et al. (2006), solar metallic- 
ity and a formation redshift of 3 to 5 as used in Grillo et 
al. (2009) are good assumptions for a galaxy of the mea- 
sured stellar mass. Thus, the IMF must not be Chabrier or 
Kroupa like unless a fraction of the de Vaucouleurs com- 
ponent is not of stellar origin, i.e. part of the dark matter 
follows the light distribution. We measure a mass to light 
ratio for the stellar component of M ^ v = (5.3l°;f ) L ^ 
using gravitational lensing and assuming a NFW-like dark 
matter halo. If we allow for a group halo at galaxy I, we 
obtain M f^ v = (6.21°'q) l ^ g b . These results again favor a 
Salpeter IMF, and are in agreement with the Fundamen- 
tal Plane results ^ = (4.8 ± 1.4) for this galaxy from 
Grillo et al. (2009). These results hold as long as the metal- 
licity is approximately solar. We compare the mass to light 
ratio, passively evolved to z=0, to Coma galaxies analysed 
in Thomas et al. (2011). We confirm their trend towards a 
Salpeter IMF, again disfavouring a Kroupa-IMF. This trend 
also is seen in Cappellari et al. (2012) for the most mas- 
sive galaxies. Their data indicate a more Salpeter-like IMF 
for high velocity dispersion galaxies. The dark to total mass 
fraction of SDSSJ 1430+4105 rises from the centre outwards, 
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giving a value of M d t ar t k = (0.4+q'Jq) at the Einstein radius. 
In this galaxy, we need a significant amount of dark matter 
in its projected centre to explain the observations. 
IV) We compare the 3d densities of total, dark and luminous 
dark matter with Coma galaxies analysed by Thomas et 
al. (2007), based on dynamical modelling, especially their 
Figure 5. Our galaxy has an effective radius of I0.96kpc 
and a velocity dispersion of 322 ± 32kms~ 1 . Concerning the 
effective radius it is most similar to the Coma galaxies CMP 
0144, CMP 4928 and CMP 2921, a cD galaxy, which have 
effective radii of 8.94 kpc, 14.31 kpc and 16.43 kpc (Thomas 
ct al. 2007) and effective velocity dispersions of 211.8 ± 
0.4kms _1 and 314.8 ± 2.9 kins -1 and w 400 kms -1 (Thomas 
et al. 2007, Corsini et al. 2008). Since our 3 dimensional 
matter densities are reliably known only between 3 kpc and 
about 6.5 kpc we decide to compare the matter densities at 
3 kpc. At this radius the matter density values of Thomas 
et al. (2007) are reliable as well. At the same time this 
radius is within the core radius (r c in Table 2 of Thomas 
et al. 2007) for all 3 CMP galaxies and thus the densities 
at this radius define the central dark matter densities in 
these galaxies. We read off dark matter and total densities 
of 6 x lO~ 2 M pc~ 3 and 3 x lO _1 M pc" 3 for CMP 0144, 
of 1.5 x lO _2 M pc -3 and 2 x IO^Mqpc -3 for CMP 4928 
and 1 x lO _1 M pc" 3 and 3 x lO _1 M pc" 3 for CMP 2921. 
For SDSSJ 1430+4105 these numbers are 4 x lO~ 2 M pc" 3 
and 3.5 x lO _1 M pc -3 for the dark matter and the total 
density at 3 kpc with fractional errors of about 25%. This 
means that the dark matter and total densities at 3 kpc for 
our galaxy and the 3 Coma galaxies are comparable, and 
that the ratio of dark to total matter density of about 1:10 
is consistent within the error with the ratios of 1:5 and 1:13 
for the non-central Coma galaxies. 
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APPENDIX A: ADDITIONAL STRONG 
LENSING MODELS 

To check for the robustness of the previously derived lensing 
results, we also examine some different strong lensing mod- 
els to the ones presented in Sees. 3.1 and 3.2. These models 
confirm the previous results without adding new implica- 
tions for the results, therefore we add these models in this 
appendix. 

Model la : To account for the environment, we include the 
galaxy group explicitly as a SIS profile centred at galaxy 
I in Tab. 3. We use a prior on the group Einstein radius 
of bg r0 u Pi p r i or = (3.6 ± 1.5)", as derived in Sec. 3.2. The 
results and 68 % c.l. marginalised errors of this SIE+GI 
(Model la) case are: b = (1.451 ,' , 2 .)", q = (0.81±^ol). 
9 q = (-17.4±*;g)° an d b group = (4.6^'j)", see Fig. Al 
for the derived parameter errors. This plot shows an anti- 
correlation of b and bgroup- This is expected since the total 
convergence needed at the position of the main lens can ei- 
ther be provided by the main lens or by the mass associated 
with GI. We also include the environment as external shear, 
as calculated in Sec Bl, see Model lb in Appendix A. This 
has only small effects on the derived parameter values. 

Model lb : Model lb is for a SIE with external shear 7, 
hence it has 1 more free parameter relative to Model la. The 
external shear priors are based on the environment models 
derived in Sec. 2.2 and Appendix B: We use 7 pr ior = (0.012± 
0.031) and 7 pr i or = (—10 ± 25)°. The marginalised errors 
are: b = (1.50+°;°!)", q = (0.8ll°;°?), 9 q = (-13.4^ 2 7 5 )°, 
7 = (0.050lo.o25) and 9 7 = (-29.6±i 5 9 8 )°. The re is a cor- 
relation present between the axis ratio q and the external 
shear 7, reflecting the fact that the shear and the ellipticity 
can compensate each other in its effects on the deflection 
angle, since both are pointing in the same direction within 
« 16°. 

Model Ila : If we add the group GI, we obtain a 
marginalised steepness value of f) = (1.71^° ^) together 
with b = (2.09i° ; 8 6 ° 9 )", q = (0.911°-°;), 6 q = (-15.3lJ° 6 2 r 
and bg r0 up = (4.51J 5)". Since there is no correlation be- 
tween ft and bg roup , the details of the environment imple- 
mentation have no systematic influence on the derived steep- 
ness of the lens mass profile. The shear and convergence pro- 
vided by the group are 7gi = kgi = 0.037 in agreement with 
our expectations. 

Model lib : A mass density profile which is flatter than 
isothermal at the Einstein radius can also be achieved by 
an isothermal mass distribution with a finite core radius. 
Therefore, Model lib is for an isothermal ellipsoid with a 
core radius (NSIE) with /? = 2 and arbitrary value for the 
core radius O c . For such a model one expects to also find a 
demagnified third image, which is not observed. We assume 
that the demagnified third image in the centre produced by 
a non-singular mass profile could be detected if its flux ex- 
ceeds 3<t of the sky+object noise in the image for the bright- 
est source pixel. We exclude a region of 0.2" in the centre 
due to residuals of the galaxy subtraction, where we have 
no limits on the image fluxes at all. We then get the follow- 
ing marginalised errors: b = (1.63lg;JJ)", q = (0.75±g;g|), 
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9 q = (-21.5+2 3) and a core radius of 6 C = (0.1ll°,;j£ 5 )". 
However, the best-fitting model is identical to Model I, i.e. 
purely isothermal. There is a linear dependency between b 
and 0c due to the definition of the non-singular profile: A 
larger core radius needs to be compensated by a larger lens- 
ing strength b to get the same enclosed mass within the 
Einstein radius. 

Model Ilia : In Model Ilia we allow for a dark matter 
component that is centred on the lensing galaxy. Here, we 
test if we can improve the modelling by combining the de 
Vaucouleurs profile with a SIE halo. For the dark matter 
SIE halo, we impose a prior on the axis ratio based on the 
Bolton et al. (2008) work of q d ark, P rior = (0.79 ± 0.12). 
For the errors, we get b = (1.431^), q d = (0.711 , , 2 ,), 
9 q , d = (-21.9+J i) and M dcV = (0.6±g;I) x lO n M . Since 
the best-fitting model has a x =11.8 and almost zero de 
Vaucouleurs mass, this implies that a de Vaucouleurs like 
mass model for the light plus a purely isothermal density 
profile for the dark matter are not compatible with the data. 
The Einstein radius of the SIE component and the de Vau- 
couleurs mass are anticorrelated, forcing the total projected 
mass within the Einstein radius to be constant. 

Model Illb : Instead of the NFW component, we can 
also model the dark matter component with a NSIE. We 
limit the core radius of this component to be between 
and 50", and use the same prior on the axis ratio as be- 
fore. For the mass of the de Vaucouleurs component, we 
get: Mdov = (8.5^0 7) x 10 n M Q . For the other parameters, 
we get (see Fig. A2): q d = (0.78±g;8|), 6 q = (-26.0l 2 j)°, 
b d = (19.2tlo:?)" and 9 C = (24.7±\H)". Again, we note 
that there is a degeneracy between b d and 8 C , emerging 
from the profile definition. Since there are no observed im- 
ages for radii larger than 2.32", this leaves the upper limit 
of the core radius 6 C totally unconstrained. Large O c make 
this dark matter distribution flat at the Einstein radius, with 
bd giving its density value. 
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Figure Al. Error estimates of the MCMC for the SIE case with 
GI (Model la), plotted are the individual points of the MCMC 
together with the 68 and 90 % confidence regions of the distri- 
bution. The crosses mark the minimum-x 2 value from Table 4. 
The bars on the axes mark the respective 68 % marginalised error 
intervals. 



Model IVa : We model the mass distribution by a de Vau- 
couleurs model with the shape parameters following the light 
profile as stated in Table 1 Therefore, the mass of the de 
Vaucouleurs component is the only free parameter in this 
model. The best-fitting light model (Model IVa) in Table A2 
has a x 2 = 118, meaning that a pure de Vaucouleurs profile 
is a bad fit to the observations. The de Vaucouleurs mass in 
this case is M doV = (15.0+. ,; 2 ,) x lO n M . This badness of 
the fit implies that there must be a (dark) mass component 
not following a de Vaucouleurs profile. 

Model Va : Here we do the same as in the Model V before: 
We combine the 2 component model (Model Illb) with an 
explicit description for the galaxy group GI. For the de Vau- 
couleurs component, we get: M de v = (11.515 5) x 1O 11 M0. 
For the other parameters, we get qd = 0.741q'iq, Q q — 
(-21.8+.^)°, b d = (5.5t 4 ;J)", 9 C = (IG.41^ 1 )", and for 
GI bg roU p = (5-lii'g)". As one can see, again there is no sig- 
nificant difference between this model's parameters and the 
one of Model Illb. As for the NFW-like dark matter halo, 
the Mdov is increased relative to Model Illb by introducing 
the group halo GI. 



Model Vb : Model Vb is motivated by the fact that for 
the preceding models (Models III and V), the axis of the 
dark matter halo is always offset from that of the light by 
about —10°, which is statistically significant on a more than 
3cr level for Models III and Illb, see Figs. 7 and A2. At 
the same time, the axis ratio of the dark matter haloes are 
consistent with the axis ratio of the stellar component, see 
Table 1. This could be mimicked by a non accounted exter- 
nal shear which is present if galaxy A is not the centre of 
the group. From the results for Models la, Ha, V and Va, 
we conclude that if we include the group explicitly as cen- 
tred on galaxy I the matter gets more aligned with the light. 
Looking at Model lb we see that using an external shear in- 
stead of GI changes the best-fitting orientation of the total 
mass distribution more towards the observed light's angle. 
So, in this Model, we combine Model Illb with the exter- 
nal shear of Model lb. In numbers, we get here: M dc v = 
(9.9+1.1) x lO n M , q d = (0.77i°;i°), 9 q = (-14.213;°)°, 
Cd = (4.ll* 7), r 200 — (I66I55)", and for the external shear 
7 = (0.038lo.o2i) and 6 7 = (-35.2l? 5 2 2 )°. We also note 
that with this improvement, the dark matter profile becomes 
more concentrated, at a level expected for galaxies. 
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Table Al. minimum-x 2 values and parameter estimates, derived with GRAVLENS for the one component isothermal and powerlaw models 
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Table A2. minimum-x 2 values and parameter estimates, derived with GRAVLENS for the two component de Vaucouleurs plus dark matter models 
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APPENDIX B: ALTERNATIVE DESCRIPTIONS 
FOR THE LENS ENVIRONMENT 

In this Appendix, we discuss 2 alternative scenarios for the 
environment, firstly a scenario in which the group is only 
consisting of its members ( "clumpy group" ) without a ref- 
erence to a group halo, secondly a scenario where the group 
is a typical group with 12 members, but centred on galaxy 
A instead of galaxy I. 



Bl Clumpy group 

A clumpy group model is obtained if all group mass is con- 
sidered to be associated to the group galaxies. We describe 
the galaxies as singular isothermal spheres (SIS) without 
truncation of their mass profiles and obtain at the position 
of the lens: 



the Faber- Jackson relation (Faber & Jackson 1976), 



clumpy group 



= ^ KSIS.n 



^clumpy group _ ^ ^ 



(Bl) 



7sis,] 



In this model, the shear and surface density at the lo- 
cation of the lens depends on the 2-dimensional galaxy dis- 
tribution and not at all on the centre of mass of the group. 
The galaxies themselves are parametrized only by their po- 
sitions and velocity dispersions <r n . The value of the veloc- 
ity dispersion <ta for SDSS J1430+4105 is taken from the 
central velocity dispersion measured by the SDSS. The es- 
timates cr n for the neighbouring galaxies are obtained from 



where i is the SDSS i band flux of the neighbours and Ia the 
flux of the lens galaxy A. The shear 7 n (d n ) and convergence 
«(d n ) for a SIS at a projected angular distance d n from its 
centre are 

c 2 d n \ D s 



^(d n ) = 7n(rfn) = 2 



with c denoting the speed of light and D^ B and D s mark the 
angular diameter distances from the lens to the source and 
from the observer to the source, respectively. The proper 
(vector) addition of these convergence and shear values 
yields a prediction of 

clumpy group _ q qj^ 

, ' ' (B2) 

^clumpy group =0 023 

The angle of the shear is —10° in the local coordinate system. 
The fact that we model the galaxies as SISs, ignoring the 
finite halo sizes which would keep the mass associated to 
galaxies limited, is not relevant, since finite halo sizes can 
only lead to lower estimates for the convergence and shear 
at the position of galaxy A. Therefore, we get an upper limit 
of the clumpy group estimates using this assumption. As we 
see in Sec. 3.2, the parameters of the lensing galaxy only 
mildly depend on the assumptions about the group as long 
as it is centred on galaxy I. To calculate the mass of this 
clumpy group, we first need the r2oo- We adopt the definition 
of Koester et al. (2007) of r2oo as a function of number 
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Figure A2. Error estimates of the MCMC for the de Vaucoulcurs+NSIE model (Model Illb), plotted are the individual points of the 
MCMC together with the 68 and 90 % confidence regions of the distribution. The crosses mark the minimum-x 2 value from Table 5. 
The bars on the axes mark the respective 68 % marginalised error intervals. The individual points of the MCMC are omitted for clarity. 



of group members. We sum up all mass contributions of 
member galaxies in Table 3 within this r2oo = 3.8' centred 
on A or I, respectively and calculate the total projected mass 
of the group within its r2oo- M200 = 5.5 x 10 14 Mq. This value 
gives an upper limit of the mass associated with the group, 
since SIS profiles for its members overestimate the densities 
of each member at large radii. 

B2 Smooth group mass distribution centred at 
galaxy A 

In principle, the assumption of the group being located at A 
could already be in conflict with the lensing observables. The 
most secure strong lensing estimate is the observed critical 
mass 7rR| in S crit = 5.43t°;i| x lO n M within the Einstein 



radius, obtained from all models in Sec. 4 consistently. We 
now can model the group - located at A as an NFW or SIS 
(see Sec. 3.1 for details) profile and estimate its projected 
mass within the observed Einstein radius. If this halo mass 
estimate exceeds the observed critical mass, the assumption 
of this group being a typical group with 12 members and 
with galaxy A as its centre is already in conflict with the 
lensing observables. 

In Fig. Bl we show the c-r2oo diagram for a NFW profile. 
The levels of grey indicate the virial M200 mass of a group 
with parameter values c and r2oo • The thick solid line marks 
the transition where the NFW group halo mass within the 
observed Einstein radius alone (without baryons and dark 
matter of the galaxy A) exceeds the critical mass, predict- 
ing a bigger than the observed Einstein radius. Therefore all 
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Figure A3. The same as Fig. 13, only for Models Illb and Va 



groups that lie above this line would - from its group halo 
mass alone - overpredict the observed total projected mass 
within the Einstein radius and cannot be centred at galaxy 
A. 

In reality, some of the observed mass within the Einstein 
radius has to be contributed by the stars, giving an even 
smaller upper limit for the dark matter mass within the 
Einstein radius. Hence, we plot the analogous curves for the 
case where the dark matter makes up only a fraction of the 
total critical mass within the Einstein radius. The dark to 
total matter fractions shown also in Fig. Bl as dash-dotted 
lines are fdark = 0.55, 0.62, 0.74. To obtain these numbers, 
we subtract the stellar mass measurements within the Ein- 
stein radius done in Grillo et al. (2009) from the derived 



lensing mass within the Einstein radius in this work. If we 
attribute the missing mass to the group dark matter halo, we 
get again upper limits for the possible group halo mass con- 
tribution within the Einstein radius, allowing us to exclude 
all groups that would exceed this upper mass limits. Grillo 
et al. (2009) fit composite stellar population models to the 
SDSS photometry of this galaxy to derive its stellar mass 
within the observed Einstein radius. We use the Salpeter 
IMF stellar masses of Grillo et al. (2009), since these give 
the highest mass in stars. Now we plot the model group with 
richness 12 in Fig. Bl to see where it resides. From Johnston 
et al. (2007), we obtain c-r2oo values of 4.22 and 848 kpc for 
a richness 12 group. Since this group therefore does not fall 
into the excluded regions of Fig. Bl, we cannot exclude A as 
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Figure Bl. This figure shows the concentration c - r2oo prop- 
erties for a NFW halo profile. The levels of grey show the virial 
masses of the dark matter haloes. Overplotted are several differ- 
ent lines: The dashed lines are the Bullock ct al. (2001) c-r200 
relation with its l<r error. This marks the area typically populated 
by galaxy groups. Further we overplot the c-r2uo values for a typ- 
ical richness N=12 group halo as found by Johnston et al. (2007) 
with its crrorbars. This shows where we expect the group halo to 
lie approximately in this plane. The dash-dotted lines mark the 
transition above which more than 55, 62 and 74 % of the observed 
critical mass within the observed Einstein radius would be made 
up by the dark matter halo of the group. All group haloes above 
this dash-dotted line in this c-r2oo plane overpredict the observed 
total mass within the Einstein radius, therefore this lines mark 
regions with excluded group haloes. Since the typical Johnston 
ct al. (2007) group halo lies below this lines, the observed critical 
mass within the Einstein radius does not exclude A as the group 
centre. The thick, solid 1.00 line marks the transition where the 
dark matter group halo alone would provide the observed critical 
mass within the Einstein radius. Hence along this line no baryons 
(or dark matter) in the lensing galaxy A would be required at all. 

the group centre from the lensing observables. This conclu- 
sion also holds in the picture where the group is modelled 
as SIS. If the group follows a SIS matter profile it has an 
Einstein radius of Oe = 3.6 ± 1.5", see Sec. 3.2. This is 
consistent within the errors with the value derived from the 
strong lensing models in Sec. 4. Therefore a typical group 
with richness 12, as seen in the vicinity of SDSSJ 1430+4105, 
does not violate the observed critical mass within the Ein- 
stein radius, nor the Einstein radius itself. Hence, galaxy A 
could also be the group centre without violating the lensing 
observables for a typical group of this richness. 



